Published online 25 January 2011 



Nucleic Acids Research, 2011, Vol. 39, No. 10 4265-4274 

doi:10.1093/nar/gkql362 



Induced topological changes in DNA 
complexes: influence of DNA sequences 
and small molecule structures 

Rebecca A. Hunt, Manoj Munde, Arvind Kumar, Mohamed A. Ismail, 
Abdelbasset A. Farahat, Reem K. Arafa, Martial Say, Adalgisa Batista-Parra, 
Denise Tevis, David W. Boykin and W. David Wilson* 

Department of Chemistry, Georgia State University, Atianta, GA 30302, USA 

Received November 5, 2010; Revised and Accepted December 30, 2010 



ABSTRACT 

Heterocyclic diamidines are compounds with anti- 
parasitic properties that target the minor groove of 
kinetoplast DNA. The mechanism of action of these 
compounds is unknown, but topological changes 
to DNA structures are likely to be involved. In this 
study, we have developed a polyacrylamide gel 
electrophoresis-based screening method to deter- 
mine topological effects of heterocyclic diamidines 
on four minor groove target sequences: AAAAA, TTT 
AA, AAATT and ATATA. The AAAAA and AAATT 
sequences have the largest intrinsic bend, 
whereas the TTTAA and ATATA sequences are 
relatively straight. The changes caused by binding 
of the compounds are sequence dependent, but 
generally the topological effects on AAAAA and AA 
ATT are similar as are the effects on TTTAA 
and ATATA. A total of 13 compounds with a var- 
iety of structural differences were evaluated for 
topological changes to DNA. All compounds 
decrease the mobility of the ATATA sequence 
that is consistent with decreased minor groove 
width and bending of the relatively straight DNA 
into the minor groove. Similar, but generally 
smaller, effects are seen with TTTAA. The intrinsic- 
ally bent AAAAA and AAATT sequences, which 
have more narrow minor grooves, have smaller 
mobility changes on binding that are consistent 
with increased or decreased bending depending 
on compound structure. 



INTRODUCTION 

Heterocyclic diamidines, and related compounds, with an 
appropriate structure can bind strongly to the DNA minor 
groove and exhibit important biological effects including 
chnical antiparasitic activity (1^). Recent results show 
that the compounds can compete with transcription 
factor binding in both the minor and major grooves 
(5,6). The compounds can also target AT-rich sequences 
in mitochondrial kinetoplast minicircle DNA (kDNA) in 
eukaryotic parasites that cause human disease (1-4), and 
they have entered Phase III clinical trials against Human 
African Sleeping Sickness (4). The kinetoplast is 
composed of maxicircles that contain the mitochondrial 
genomic DNA and thousands of 1000-2000 bp catenated 
DNA minicircles that code for RNAs that guide editing of 
the maxicircle RNA transcripts (7). The kDNA minicircles 
have a number of unusual properties and structural 
features, such as short, curved A-tract sequences that are 
in phase with the helix repeat so that they always appear 
on the same side of the double helix and cause global 
bending of the helix (7-9). Given these unusual structural 
features, it seems hkely that topological changes in the 
DNA minicircles induced by minor groove binding com- 
pounds can lead to kinetoplast destruction during opening 
of minicircles for transcription and replication and, thus, 
play a part in the antiparasitic activity of the compounds 
(1). Binding-induced changes in DNA topology are also 
important in the regulation of transcription of many 
genes, and modulation of bending can enhance or 
reduce transcription (10-12). Small molecules that can 
affect DNA bending in different ways at selected sites 
could, thus, provide tunable controls on the expression 
of the target genes. 
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In order to test the hypothesis that compound- 
induced topological changes in kDNA are important for 
their antiparasitic activity, a polyacrylamide gel electro- 
phoresis (PAGE) ligation ladder assay (13-17), with 
phased AAAAA bent sequence and ATATA relatively 
unbent sequences, was established to determine whether 
heterocyclic compounds can cause topological changes 
in DNA (18). Out of phase hgation ladders of the same 
AT sequences were used as controls. Ligation ladders with 
repeated identical binding sites in phase with the helix 
twist ampHfy and simplify the analysis of compound topo- 
logical effects at a specific sequence and different binding 
sites can be evaluated in the same context. The results 
in the initial assay showed quite clearly that the A-tract 
sequences were shghtly straightened by DB75, while the 
compound significantly bent the alternating AT DNA 
sequence (18). Compounds with single atom changes in 
DB75 were also evaluated by replacing the furan with 
thiophene (DB351) and selenophene (DB1213) groups 
and compound-specific effects were observed. In these 
experiments the gel electrophoresis was done with the 
compound included in the gel at a concentration to 
maintain saturation of binding sites during the experiment 
(Method 1). The method of putting compound in the 
gel to help prevent dissociation has been successfully 
used with other compounds (19,20). The heterocyclic 
diamidines investigated here are positively charged at the 
typical buffer pH of 8.3, and thus migrate in a direction 



opposite that of the DNA during electrophoresis. 
Including the compound in the gel (Method 1), however, 
hmits the experiment to analysis of one compound at one 
concentration per gel. 

The interesting DNA topological changes observed on 
binding of DB75 and analogs suggested that it would be 
worthwhile to evaluate additional compounds as well as 
DNA sequences to determine the generahty of the obser- 
vations. Such studies will also help to estabhsh relation- 
ships for compound structural effects on the topology of 
kDNA and other biologically important sequences. In 
order to increase the number of compounds and DNA 
sequences that can be evaluated per gel, a method that 
does not use compound in the gel is needed. Here a dif- 
ferent PAGE method (Method 2), without compound in 
the gel, has been tested. Method 2 was tested with com- 
pounds of varied structure (Figure 1), and binding con- 
stants were determined to estabhsh the limits for 
application of the method. The PAGE method without 
compound in the gel is able to qualitatively reproduce 
the more quantitative results with Method 1, and with 
Method 2 it is possible to study several compounds at 
different ratios to DNA sites with different DNAs on 
one gel. 

Compound selection 

Compounds were selected for testing in the screening 
assay based on systematic modifications to each of the 




Figure 1. Compounds grouped by structural similarity. 
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molecular units of DB75 (Figure 1). The compounds are 
grouped based on common structural motifs (Figure 1). 
The central furan of DB75 was replaced with a pyrrole 
(DB262) and from that structure the terminal amidines 
were changed to reversed amidines (DB884). One phenyl 
was changed to a benzimidazole (DB293) and from that 
compound the furan was changed to an imidazole 
(DB940), or the other flanking phenyl was replaced with 
a furan (DB1003). Another group had one amidine 
replaced with an amide (DB1932 and DB1933) such that 
they were monocations, not dications. DB185 is a different 
structure that is related to the well-known bisbenzimi- 
dazole, Hoechst 33258 (Figure 1). A final group of 
compounds has one phenyl of DB75 changed to 
benzimidazole and the furan converted to a phenyl. In 
this group DB911 has the substituents in a meta position 
at the central phenyl to give a classically curved minor 
groove binding shape, while an isomer has the central 
phenyl para substituted (DB921) to give a relatively 
linear compound. Two analogs of DB921, one with the 
addition of a hydroxyl group (DB985) and one with the 
position of the terminal phenyl amidine group in a meta 
position (DB1055), are included in this group. 

MATERIALS AND METHODS 

DNA binding compounds 

Compounds in Figure 1 were synthesized by previously 
described methods (21-26), while the synthesis of 
DB940, DB1003, DB1932 and DB1933 will be pubhshed 
elsewhere. All the compounds were prepared in a 1-mM 
stock solution and stored at 4°C away from hght. 



DNA oligonucleotides 

Two 21-nier duplexes (AAAAA and ATATA in Figure 2) 
from Tevis et al. (18) and two new duplexes (TTTAA and 
AAATT in Figure 2) were used in this study. The marker. 
Ml, is also a 21 mer with the same base composition 
as the test duplexes, however; there is no AT binding 
site (18). Marker, M2, is a commercially available 20 bp 
ladder (Bayou Biolabs, Harahan, LA, USA) with a double 
intensity band at 100 bp. High-performance liquid chro- 
matography (HPLC) purified and lyophilized DNA oligo- 
mers were purchased from Integrated DNA Technologies, 
Inc. (IDT, Coraville, lA, USA). The concentration of 
single-strand DNA was brought to ~ 1 mM using water. 
The exact concentration was then determined using the 
extinction coefficient calculated by the nearest neighbor 
approximation. Then the complementary strands 
(100 |iM) were combined in a 1:1 ratio based on the 
determined concentrations and annealed in Ix ligation 
buffer (New England BioLabs, Ipswich, MA, USA) con- 
taining 50niM Tris-HCl, lOmM MgClz, lOmM 
dithiothreitol, 1 niM ATP and 25 mg/ml bovine serum 
albumin. 

Ligation ladders 

The ligation ladder protocol is the same as reported by 
Maher et al. (15,27,28) and Tevis el al. (18). Annealed 
duplexes (100 |iM) were 5'-phosphorylated using T4 poly- 
nucleotide kinase (New England Biolabs) and then ligated 
using T4 DNA hgase (New England Biolabs) in the buffer 
provided. The kinase reactions were 200 jil of 6 jiM DNA 
with 80 U of T4 polynucleotide kinase. The reaction 
was incubated for 30min at 37°C followed by enzyme 



1 AAAAA 

5'- GAAAAACGCCGCAAAAACGCC 



GGCTmTGCGGCGTTl 
2 TTTAA 



ITGC -5' 





A5 


TA 




ATATA 




M2 M) Mi4^ 


0 2:1 


0 2:1 


0 2:1 


0 2:1 





5'- GTTTAACGCCGCTTTAACGCC 



GGCAAATTCCGGCGAAATTGC -5' 
3 AAATT 

5'- GAAATTCGCCGCAAATTCGCC 



GGCTTTAAGCGGCGTTTAAGC -5' 
4 ATATA 

5'- GATA1ACGCCGCA1ATACGCC 




GGCTATATGCGGCGTATATGC -5' 
5 Ml 

5'- GACAGACACAAGACACAGACC 
GGCTGTC TGTGTT OGTGTCT-S' 

Figure 2. Sequence of DNA duplexes used (left) and 8% PAGE gel of DB75 with the duplexes (right). The colored arrows indicate the 189-bp 
duplex for each ladder. The red arrow of M2 indicates 100 bp, while the red arrow of M] indicates 189 bp of the random sequence marker. 
Electrophoresis was done for 2h 20min at 200 V and compound was not included in the gel. For each duplex, 0 indicates DNA that is not 
incubated with compound prior to electrophoresis, and the ratio 2:1 (compound to binding site) indicates DNA incubated with compound prior 
to electrophoresis. The brightest bands near the top of the lanes are circular products. 
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deactivation for 30min at 65°C. Ligation reactions were 
200 (il using 1200 U of ligase and ligation was done at 
room temperature for 20min followed by enzyme deacti- 
vation for 25min at 65°C. 

Gel electrophoresis 

Ligation ladders were separated using 8% native poly- 
acrylamide (1.5 mm thick, 20 cm long) prepared from a 
40% acrylaniide stock solution (29:1 bis-acrylaniide:- 
acyrlamide; HMD, Gibbstwon, NJ, USA) in Ix TBE 
buffer (0.089 M Tris, 0.089 M boric acid, 2.0 niM EDTA, 
pH 8.3). Each sample had a final volume of 20 |.il, which 
contained the following: 10 ^l of 2^M DNA, 4^1 of 
6x load dye (Promega, Madison, WL USA), the volume 
of compound needed for the desired ratio of compound 
to binding site and double distilled water to bring to the 
final volume, if needed. For the 1:1 and 2:1 ratios 3 and 
6^1 of compound (13.3 |iM) were added, respectively. 
Electrophoresis was done at 200 V (lOV/cm) for 140 min 
at 25°C ± 0.1°C in Ix TBE buffer. The quite different 
observed mobihties with compounds in adjacent lanes 
indicate that when using the 140-min time period, there 
is no significant interference from inter-lane diffusion 
when different compounds are used in the gel lanes. The 
lack of interference is also supported by the observation 
that compounds have the same mobihties when studied 
alone (Figures 2 and 3) or on a gel with other compounds 
(Figure 4). 

After electrophoresis, gels were stained with SYBR 
Gold Nucleic Acid stain (Invitrogen, Carlsbad, CA, 
USA) at the concentration and time recommended by 
the manufacturer. The gels were imaged using an 
UltraLum Omega 10 gD Molecular Imaging System 
(UltraLum, Claremont, CA, USA). 



Calculating /?l 

The analysis of topology changes was done with 
ImageQuantTL (Amersham Biosciences, Piscataway, NJ, 
USA). From the base pair (bp) assignments the parameter 
i?L, where i?L = -f-apparent/^actuab was Calculated 
(18,27,28). The mobility of each band of the standard. 
Ml, was determined and correlated with its length, in bp 
(Lact) and molecular weight. The mobihties of each band 
of the test DNAs and their complexes with diamidines 
were next determined and correlated to apparent bp 
lengths (^apparent) ^nd molccular weights based on the 
standard values. The i?L value for each unbound DNA 
duplex was normalized to 1.00 so that the changes 
induced by compound binding could be more easily 
compared within each sequence. 

Biosensor-surface plasmon resonance studies 

AAAAA, ATATAT, AAATT and TTTAA hairpin 
duplexes with 5'-biotin groups were purchased from IDT 
(Integrated DNA Technology) with HPLC purification 
and mass spectrometry characterization. The oligomer 
sequences are all as follows: GCGWWWWWCG TCTC 
CGWWWWWCGC, where the W sequences are AAAA 
ATTTTT, ATATAT ATATAT, AAATT AATTT, TTT 
AA TTAAA and the hairpin loop is underhned. A 1-mM 
stock solution of each compound was prepared in water 
and diluted to working concentrations with buffer imme- 
diately before use. Surface plasmon resonance (SPR) 
measurements were performed with a four-channel 
Biacore 2000 optical biosensor system by using pubhshed 
procedures (29,30). For the duplex binding deter- 
minations, the 5'-biotin-labeled DNA samples were 
immobilized in separate flow cells on streptavidin-coated 
sensor chips (Biacore SA) with one flow cell left blank as a 




Figure 3. (A) PAGE (8%) of c/s-ligated sequences with DB921 not in the gel. Electrophoresis was done as in Figure 2. For each duplex, 0 indicates 
DNA that is not incubated with compound prior to electrophoresis, and the ratios (coinpound to binding site) indicate DNA is incubated with 
compound prior to electrophoresis. The colored arrows indicate 189 bp for each duplex. The red arrow of M2 indicates 100 bp, while the red arrow of 
Ml indicates 105 bp of the random sequence marker. (B) PAGE (8%) was done for 2h 20 min at 200 V and compound was included in the gel. 
Experiments were done with the same cis sequences as in (A) and with traii.s-ligated DNA sequences of AAAAA (tAs) and ATATA (tAT). The 
colored arrows indicate 189 bp and the double intensity band of M2 is 100 bp. 
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Figure 4. PAGE (8%) of the DNA sequences (Figure 2) with multiple DB compounds. (A) The duplexes AAAAA and ATATA with five selected 
compounds, and (B) TTTAA and AAATT with the same five compounds. For each duplex, DNAs are incubated with compound prior to electro- 
phoresis at a ratio of 2:1 (compound to binding site). The colored arrows indicate 189 bp for each duplex. The red arrow M] indicates 105bp of the 
random sequence marker. Electrophoresis was done as in Figure 2 and compound was not included in the gel. 



control. The SPR experiments were performed at 25°C in 
filtered, degassed lOmM cacodylic acid buffer (pH 6.25) 
containing 100 mM NaCl, 1 mM EDTA and 0.005% sur- 
factant P20. Compound solutions were prepared by serial 
dilution from the stock solution and injected from 7-mm 
plastic vials with pierceable plastic crimp caps at the flow 
rate of 25 ^il/min. The compound was then replaced by 
buffer flow to monitor dissociation of the complex. The 
reference response from the blank cell was subtracted 
from the response in each cell containing DNA to give a 
signal (response units, RU) that is directly proportional to 
the amount of bound compound. Sensorgrams, RU versus 
time, at different concentrations for binding of each 
compound to DNA were obtained, and the RU in the 
steady-state region was determined by linear averaging 
over a selected time span. The number of binding sites 
and the equilibrium constant were obtained from fitting 
plots of RU versus Ctvee- Binding results from the SPR 
experiments were fit with either a single site model 
{Kj = 0) or with a two-site model: 

r = {K\ X Crree+2 X Ki X K2X Cfj.ee)/ 
(l+KiX Cfree+2 xKixKjx Clj 

where ;• represents the moles of bound compound per mole 
of DNA hairpin duplex, Ki and K2 are macroscopic 
binding constants and Cfree is the free compound concen- 
tration in equihbrium with the complex. 

RESULTS 

Topological screening methods: comparison of two assays 

In Figure 2 the A5 sequences show increased mobility 
(green arrows), while the ATATA sequences have a 
reduced mobility (blue arrows) relative to the same 



DNAs without DB75, and the results are very similar to 
those previously reported for the same sequences with 
the compound in the gel (18). Two new duplexes, TTTAA 
and AAATT, are also included in Figure 2 and the 
pink arrows show that DB75 does not have a detectable 
topological effect on TTTAA under these conditions, 
while the orange arrows show that DB75 increases the 
mobihty of AAATT to an even greater extent than with 
A5. There are two alternative explanations for the 
mobihty changes observed for DB75 and the other com- 
pounds in this study; the compounds could be changing 
the bend angles of the DNAs or they may be slightly 
changing the helical repeat, resulting in improved or 
decreased phasing of the AT repeats. Because the trans 
A5 sequence is already significantly dephased, an 
increase in the mobility of the A5 ladder is probably due 
to a decrease in the curvature of the A5 A-tracts. In the 
same manner the dramatic reduction in mobihty of the 
ATATA sequence is highly unlikely to result from a 
phasing change because the sequence is straight no 
matter whether it is phased or not (18 and cis-trans in 
Figure 3). It seems hkely then that bending rather than 
dephasing of helical repeats is responsible for the observed 
changes in mobility for these DNAs and compounds. 
Additional experiments with 20- and 22-bp duplexes 
would be required in order to definitively differentiate 
between phasing and bending or straightening of the 
AT sequences to account for the mobihty changes. 

The brightest bands near the top of each lane are 
circular products that appear for all DNAs as the molecu- 
lar weight increases. The circular bands appear at lower 
molecular weights for the bent DNAs, AAAAA and AAA 
TT than for the more hnear duplexes. As can be seen in 
the figure, the compounds have very little effect on the 
mobihties of the circular DNAs. As compounds are 
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added, the linear bands can change markedly in mobihty 
while there are no detectable changes for the circular 
DNAs. The potential significance of this for biological 
activity is presented in the 'Discussion' section. 

Figure 3 shows an additional comparison of gels done 
with DB921 by Method 1, which has DB921 in the gel 
(Figure 3B), and the new Method 2, where DB921 is 
added only to the separate DNAs and is not in the gel 
(Figure 3A). DB921 has a very different structure from 
classical minor groove binders with a more hnear shape 
that incorporates an interfacial water molecule to bind in 
the minor groove (31). As can be seen, DB921 reduced the 
mobihty of aU four of the c«-hgated duplexes of interest 
by similar amounts in both assays and the results are quite 
different from those with DB75 (Figure 1; 18). This 
example shows that both methods give equivalent results 
for compounds of quite different structure. Interestingly, 
for the TTTAA and ATATA sequences there is a 
smearing of the bands in Figure 3A that is not present 
in Figure 3B and this smearing appears to be the result 
of partial dissociation of the DNA-compound complex 
during the gel experiment. In order to help clarify this 
difference, binding constants (see below) were determined 
for DB921 and show that it binds more weakly to TTTAA 
and ATATA than to AAAAA and AAATT. These results 
explain the smearing of DB921 bands in Figure 3 A 
because it will more easily dissociate from the TTTAA 
and ATATA DNAs during the gel experiment. Note, 
however, that the general effects of DB921 on all the 
DNAs are comparable in both methods and can be 
determined from the screening method. 

For comparison with our previous results (18), Figure 3 
also includes the trans AAAAA and ATATA sequences. 
With DB75 we observed much smaller effects on both the 
bent AAAAA sequence and the straight ATATA 
sequence. DB921 reduces the mobihty of the trans 
duplexes, but the effects are much less than with the cis 
sequences, as expected. 

Effects of compound structure on DNA topology 

In Figure 4, PAGE results for two DNAs and multiple 
compounds at a fixed ratio are compared. Topological 
effects of the compounds used in these gels agree with 
the results in Figures 2 and 3 and Supplementary Data. 
The screening method with multiple compounds, thus, 
allows for rapid testing and direct comparisons of topo- 
logical changes for compounds in a hbrary. In Figure 4 the 
2:1 ratio is used, because smaller effects are seen at a 1:1 
ratio and going to a 3:1 ratio did not generally cause large 
changes in mobility with these DNA sequences and com- 
pounds. The PAGE results for all the other compounds 
in Figure 1 can be seen in Supplementary Data. In the gels 
for DB1003 and DB185 (Supplementary Figures SM4 and 
SMIO, respectively) additional ratios were included and 
none of the extra ratios changed the conclusions. To 
accommodate the extra ratios, the ATATA sequence 
was omitted for DB1003 and DB185. These two com- 
pounds were selected for this test because DB1003 is a 
relatively weak binder (32), while DB185 is a very strong 
minor groove binding agent (33). 



Using /Jl to compare topological changes 

i?L values for each DNA and complex were calculated 
relative to the mobility of the standard as described in 
the 'Materials and Methods' section. Table 1 compares 
i?L values of several compounds that were determined in 
both assays. The topological changes shown in Table 1 are 
similar by both assays, however; for some compounds and 
sequences the effects on i?L are larger by Method 1. This 
difference is not unexpected from the results in Figure 3 
and as noted above, is due to partial dissociation of 
the complex in the screening assay (see binding results 
below). The trend and general magnitude of the changes, 
however, can be determined without compound in the gel. 
For quantitative determination of bending angle changes 
on compound binding, however. Method 1 should be 
used. 

To assist in comparing results for different DNAs, all 
i?L values were normalized by dividing them by the i?L 
value of the free DNA, and the results are plotted in 
Figure 5. A compound with an Rj^ value greater than I 
indicates the sequence had lower mobility than the free 
DNAs and is consistent with bending, while a value 
below one is consistent with straightening of the 
sequence. Interestingly, the ATATA sequence was bent 
by all compounds tested although the magnitude of 
bending varied with compound structures. DB921 and 



Table 1. Ri^ values for both assays for each DNA sequence 



AAAAA TTTAA AAATT ATATA 



Drug 


C 


No C 


C 


No C 


C 


No C 


C 


No C 


DB75 


0.77 


0.82 


1.20 


1.01 


0.75 


0.86 


1.31 


1.19 


DB262 


0.92 


0.95 


1.25 


1.06 


0.76 


0.81 


1.48 


1.27 


DB884 


1.22 


1.11 


1.16 


1.08 


1.05 


1.03 


1.19 


1.13 


DB921 


1.1 


1.12 


2.1 


1.33 


1.09 


1.08 


2.5 


1.9 



C indicates compound in the gel matrix (Method 1), and No C indicates 
no compound in the gel matrix (Method 2). The Method 2 comparison 
uses the values at the 2:1 ratio of compound to binding site. 




Figure 5. Histogram showing values for all compounds with 
AAAAA, ATATA, TTTAA and AAATT sequences. The blue repre- 
sents AAAAA, red ATATA, the green TTTAA and the purple 
AAATT. Each group represents a structural group from Figure 1. 
All the data shown are for a 2:1 ratio of compound to binding site. 
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DB911 show the largest induced bending, followed by 
DB262 and DB75. All three of the monocations caused 
a small reduction of mobility and a small bend with the A 
TATA sequence. 

The changes in i?L for the AAAAA sequence were both 
positive and negative and were generally smaller than 
the changes with ATATA. The following compounds 
cause a decrease in mobihty, consistent with increased 
bending of the AAAAA sequence: DB1933, DB921, 
DB1003, DB911 and DB884. DB75 and DB185, on the 
other hand, show an increase in mobihty and DB262, the 
pyrrole analog of DB75 increase mobihty but to a lesser 
extent. Interestingly, the monocation, DB1933, causes the 
largest decrease in mobihty in AAAAA. DB884, with a 
reverse amidine cation, also decreased the mobility of the 
intrinsically bent AAAAA sequence. The compounds 
DB293, DB940, DB985, DB1055 and DB1932 did not 
have a significant topological effect on the AAAAA 
sequence. In general, the results with AAATT are 
similar to those with AAAAA. 

For the TTTAA sequence the group containing the 
phenyl-benziniidazole modification (Figure 1) shows the 
largest effect on bending. DB921 and DB911 increase 
bending more than DB1055 that bends more than 
DB985. DB262 increased the degree of bending more 
than DB75. For the benzimidazole modification group, 
DB293 caused the largest bending effect with TTTAA 
and the Aw-benzimidazole DB185 shows a similar effect. 
The effects with DB75, DB940, DB1003 and DB1932 are 
not significant for the sequence. 

Biosensor-SPR determination of equilibrium binding 
constants {K) 

Binding assays were conducted for the compounds from 
the topological studies to answer two important questions: 
(i) what is the minimum K value that is necessary before 
Method 2 can be used without compound in the gel and 
(h) is there any correlation of the K and i?L values. 
Because three biotin-labeled DNAs can be immobilized 
on four-channel Biacore streptavidin-coated sensor chips 
with one flow cell left as a blank (29,30), two sensor chips 
were used to evaluate K values for the four DNA 
AT-binding sequences used in the topology-ligation 
assays. This provided repeat values for two DNAs to 
give a better evaluation of reproducibility of K values 
determined with completely different experimental sensor 
chips, samples, etc. The variation in K is <10% in the 
repeat assays. 

Example sensorgrams for DB262, which has similar 
structure and DNA binding to DB75, show that the 
compound has relatively fast association and dissociation 
kinetics for binding to ATATA and AAATT (Figure 6A 
and B). In this case K values can be determined by fitting 
the kinetic constants for association (/Ca) and dissociation 
{k^ or from the steady-state response where the rates of 
association and dissociation are equal (Figure 6A and B). 
The steady-state method is not subject to any mass trans- 
port issues at the sensor surface, and steady-state RU 
values for DB262 and three DNA sequences are plotted 
in Figure 6C. values were determined by fitting the plots 
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Figure 6. Example biosensor-SPR sensorgrams (instrument response in 
RU values versus time) for binding of DB262 to (A) ATATA and 
(B) AAATT at 25°C. The unbound compound concentrations in the 
flow solutions range from 0.005 nM in the lowest curve to 1 |xM in the 
top curve. Experiments were conducted in cacodylic acid buffer with 
0.1 M NaCl. (C) SPR binding plots for DB262 in the presence of 
ATATA (top), AAATT (middle) and TTTAA (lower). The RU 
values from the steady-state region of SPR sensorgrams are plotted 
versus the unbound compound concentration. The RUs and the con- 
centrations of DB262 are from repeat experiments. The lines in the 
Figure are for one or two-site fitting ('Methods' section). K values 
for AAATT were obtained by using a one-site model and for 
ATATA and TTTAA with a two-site model to get A", and K2. K2 
was much lower than K\ and only K\ values were used in the compari- 
son plot in Figure 7. 

in Figure 6C to one or two site functions as described in 
the 'Materials and Methods' section. In ah cases the Kx 
values are significantly larger than the K2 values when a 
two-site fit was required, for example, with the ATATA 
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and TTTAA DNAs (Figure 6C). A simple two-site 
binding function can be used in these cases because the 
primary binding is for complex formation at the AT target 
sequences while the much weaker, secondary binding is 
generally a non-specific electrostatic interaction (34). 
Results for DB921, which has slower association and dis- 
sociation rates, are shown in Supplementary Figure SMI 5, 
Supplementary Data. K values for DB921 and related 
compounds were determined by global fitting of the 
sensorgrams to determine k^^ and values as a function 
of the compound concentration in the flow solution. 
DB921, DB985 and DB1055 are generally quite strong 
minor groove binding compounds. For visualization, 
primary binding constants (^i) are compared for repre- 
sentative compounds as a histogram in Figure 7, and the 
actual K values are tabulated in Supplementary Figure 
SM15. 

As can be seen, there is a large variation in binding 
constants for the compounds with different DNAs. In 
general, the highest binding constants were observed 
with the AAATT site and the lowest with the reversed 
sequence TTTAA. Significant outliers in this pattern are 
DB1003 and DB1932, which bind best to TTTAA as 
expected for compounds with that structure (32). The 
compounds with the largest K values in this group are 
DB921 and analogues, but they also have relatively low 
K values for TTTAA and the highest values for AAAAA 
and AAATT. The two pyrrole compounds, DB262 and 
DB884, have quite large values for AAATT relative to 
their ^s for the other three DNAs. Most compounds 
have larger K values for the non-alternating, AAAAA, 
than for the alternating ATATA sequence. Interestingly, 
the diamidines, DB75 and DB262, with a diphenyl- 
5-member ring structure, have larger values for ATATA. 
As will be discussed below and as shown in Supplementary 
Figure SMI 6, there is no direct correlation of binding 




Figure 7. Histograms of equilibrium constants. A' (M ') for representative 
compounds binding to ATATA, AAAAA, AAATT and TTTAA. 



constants with the mobility changes for the compounds 
and DNAs. 



DISCUSSION 

Experiments with ligation ladders from DNA sequences 
with AAAAA and ATATA binding sites have shown that 
DB75 and some close analogs induce topological changes 
in DNA that are sensitive to DNA sequence/structure and 
to compound structure and properties (18). At present 
there is considerable knowledge of the effects of protein 
binding on DNA topology but much less information 
about small molecules that target DNA. Such compounds 
are being designed for therapeutic applications that 
range from antiparasitic to anticancer activities (1-4), 
and in order to develop these agents in a rational way 
we need to understand the topological effects at the 
target sequences. Topological changes on binding of the 
heterocyclic diamidines of Figure 1 to the mitochondrial 
kinetoplast DNA of disease-causing eukaryotic parasites 
may be a key step in destruction of the kinetoplast and the 
biological mechanism of action of the compounds (3). To 
clarify this question, more topological information on the 
effects of the compounds on kinetoplast target sequences 
is required. 

Method 1, with compound in the gel to prevent the 
cationic compounds from dissociating from DNA during 
electrophoresis, has been developed for quantitative evalu- 
ation of compound-induced topological changes on 
complex formation (18). This method, however, allows 
only one compound at one concentration (ratio per 
DNA binding sites) on each gel. In this work. Method 
2, in which the compounds are added to DNA at different 
ratios but are not in the gel, was tested. This method 
allows the possibiUty of different compounds on each gel 
and on a single gel with 20 lanes, for example two ratios 
of different compounds can be tested with up to four dif- 
ferent DNAs, free DNA and appropriate standards. 
Figure 3, with one compound, and Table 1 with four com- 
pounds and four DNAs compare results by Methods I 
and 2. In Table 1 the changes observed by Method 1 are 
generally greater than those with Method 2, as expected if 
there is some compound dissociation and loss during elec- 
trophoresis. Partial dissociation is supported by results in 
Figure 3, where some smearing is seen in gel A (Method 2) 
with the TTTAA and ATATA sequences, but no smearing 
is observed in gel B (Method 1). The smearing effect 
is observed with compounds and DNAs where there is 
a large difference in Rj^ with a relatively low K value. 
The smearing, for example, occurs with the TTTAA and 
ATATA sequences, and DB884, which has lower binding 
constants with these DNAs than with AAAAA and 
AAATT (Figure 7 and Supplementary Figure SMI). The 
encouraging aspect of these comparisons, however, is that 
in all cases the direction of the changes and general 
relative magnitude is reported identically by the two 
methods. These results indicate that Method 2 can be 
used for qualitative analysis of induced topological 
changes by small cations that target DNA. For com- 
pounds that are found to be particularly interesting from 
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these results, more quantitative, but slower, studies can be 
done by Method 1 as necessary. 

With these encouraging results for Method 2, several 
groups of minor groove binding compounds of quite dif- 
ferent structures and properties (Figure 1) were selected 
for a more extensive, comparative analysis of induced 
topological effects on complex formation. The set of com- 
pounds and DNAs is satisfactory to determine if the 
induced changes are dominated by either compound 
or DNA structure and properties, and also to look for 
any patterns in the results. A qualitative conclusion 
from the results that are summarized in Figure 5 for i?L 
and Figure 7 for K is that there is no direct correlation 
between the magnitude and direction of the induced topo- 
logical changes and the binding affinity, as shown in 
Supplementary Figure SMI 6. As will be described 
below, the observed topological changes appear to be 
best correlated with DNA intrinsic minor groove width 
as weU as with compound shape and properties. 

The AAAAA minor groove is intrinsically narrow with 
ordered water and cations stabilizing the close approach 
of phosphate anions across the groove (35,36) and as a 
result, this sequence is the most intrinsically bent in our set 
of four DNAs (27,28). The groove is of the correct width 
to bind most compounds of Figure 1 with little change 
in width. The compounds can cause a shght increase in 
bending by stabilizing a slightly more narrow minor 
groove with increased bending into that groove. 
Decreased bending, which can be caused by the opposite 
effect of a slight widening of the groove, might be gener- 
ally expected, and for this compound set it is observed for 
most complexes. Effects on the AAATT sequence, which 
also has A-tract properties and significant bending for the 
free DNA, are generally similar to those for AAAAA and 
can be either in the direction of increased or decreased 
bending. 

The ATATA sequence is the least intrinsically bent 
DNA in our set of four (18,37-42) and this DNA has a 
wider groove than AAAAA or AAATT. The groove is too 
wide to bind tightly to the compounds of Figure 1 without 
some structural changes. Minor groove binding that de- 
creases the groove to better fit the general width of the 
compounds should bend the DNA into the minor groove. 
All compounds of Figure 1 whether mono or dications, 
with or without a benziniidazole group, and with amidines 
or a reversed structure (DB884) bend ATATA indicating 
that topological effects with this sequence, are dominated 
by the DNA structure. The benzimidazole-biphenyl 
isomers, DB921 and DB911, cause the most bending, 
and the related compounds, DB985 and DB1055, cause 
a smaller but significant amount of bending. The 
biphenyl motif of these compounds is intrinsically 
twisted by 30^0° and this may explain the relatively 
large bending effects of the compound family. Replacing 
the central phenyl by a furan, DB293, or with an imid- 
azole, DB940, results in a planar system with smaller 
induced bends in ATATA. The diphenyl furan, DB75 
and pyrrole, DB262, compounds are also planar, 
however; and cause significant bending. Clearly 
compound shape and perhaps other features are import- 
ant in addition to compound twist. The clear pattern with 



ATATA is that all minor groove binders studied to date 
cause bending of that sequence. The magnitude of 
bending, however, depends on the interaction of the 
compound structure and chemical groups with the ATA 
TA groove. The TTTAA sequence also has a wider groove 
than the A-tract sequences and it has a much lower degree 
of bending (41,42). Compounds also generally show 
bending on complex formation with this DNA, but the 
bending magnitude is lower than with ATATA. This 
could be due to the fact that compounds, such as 
DB1003 and DB293, have been shown to bind to the 
TTAA sequence as a dimer (32), which will maintain a 
wider groove. 

It is a very interesting result that the circular DNA 
bands in Figures 2-4 (the bright bands in the upper part 
of the gels, see Figure 2) do not change mobihty by any 
significant amount when compound is added. This is true 
even though they have the same sequence as the hnear 
DNAs and even when the mobilities of the hnear bands 
change by a large amount. Clearly constraining the ends 
of these sequences is the dominant feature of their shape 
and mobility in PAGE. Because the kinetoplast DNAs 
are circular, this observation suggests that the compounds 
would have relatively little effect on structure until the 
DNAs were opened for transcription or replication by 
mitochondrial enzymes. At this point their topology 
could change significantly and closing the DNAs to 
reform the kinetoplasts would be inhibited. Given the 
thousands of minicircles in each kinetoplast (7), even a 
small amount of inhibition could lead to very significant 
destruction of the kinetoplast structure and network. 

In summary, the major features of both induced topo- 
logical changes and of compound binding constants for 
minor groove complex formation are dominated by the 
groove characteristics of DNA AT sequences. There are, 
however, major variations within a DNA sequence for 
changes in compound structure and properties. The 
screening method described in this report can now be 
used with a quite large set of DNA binding agents to de- 
termine correlations with molecular structure and induced 
topological changes. 
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